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),  An  Evaluation  on  Damping  Property  of  Non-linear  Shimmy 
Damper  and  Equivalent  Linearization  Methods  \ 

J 

Nanchang  Aircraft  Manufacturing  Co. ,  Su  Kaixin 


Abstract 


In  £his  paper  ^t he  formula  of  requisite  linear  damping  for 
overcoming  landing-gear  shimmy  4s  der-i-ved^according  to  the 
point-contact  theory.  However,  the  damping  property  of  the 
hydraulic  shimmy  damper  is  mostly  non-linear.  Two  methods 
are  presented  to  treat  this  case:  the  tangent  method  and  the 
minimum  equivalent  linear  damping  method.  Their  correlation 
is  also  found.  It  is  demonstrated  by  computation  that  both 
methods  are  effective  and  reliable.  "  ft 

I.  Simplification  of  Theoretical  Equation  for  Single  Wheel 
Shimmy  without  Inclination 

Based  on  the  following  assumptions  that  (1)  the  support  (or 
the  directional  axis)  is  vertical  and  rigid,  (2)  a  static 

manuscript  received  on  March  4,  1987 


rigidity  model  is  used  to  treat  tire  distortion  and  there  is 

pure  rotation  and  no  skidding,  (3)  the  entire  system  is 

» 

linear,  and  (4)  it  is  moving  at  a  uniform  speed  and  the 
vertical  load  is  also  constant,  then  the  equation  of  motion 
for  shimmy  is: 


1V-—  +hvi--aXt-b<P=0 


ds‘ 

d(0  +  <P ) 
ds 


ds 

—a. :.+0$>=O 


0) 


where 

I  -  rotational  moment  of  inertia  of  the  shimmy  part  with 
respect  to  the  directional  axis, 
t  -  distance  of  stability, 

h  -  damping  moment  coefficient  centered  around  the 
directional  axis, 

V  -  speed, 

s  -  rolling  distance  along  the  direction  of  motion, 

0  -  angle  between  the  diametral  plane  after  the  wheel  yawing 
and  that  of  the  wheel  in  neutral  position, 
a,  (5  -  motion  constants  of  the  wheel, 

a,  b  -  lateral  rigidity  and  torsion  rigidity  of  the  wheel, 

A — -  variables  introduced  in  the  point-contact  theory  to 
describe  torsion  distortion  and  lateral  distortion  of  the 
tire.  Pis  the  yawing  angle  of  the  tangent  at  the  contact 
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center  and  counterclockwise  is  positive.  X  is  the  lateral 

shift  at  the  contact  center  and  it  is  positive  toward  the 

» 

right  (see  Figure  1). 


Let 


0  =  Ae’‘  1 

i 

<p=Be‘  *  1  (2) 

A  =  Ce'  *  ! 


Substitute  it  into  equation  (1)  and  let 


flt=ai>+oy St  'i 

al’ssabt~f-afiti  +  hVa  j 

at«hV)3+b+ot*  +  JaV*  !  (3) 

aj«IV*/?+ hV  | 

a*=JV* 

He  then  get  the  following  eigen  function: 

a4*4  +  ajx,  +  fl1**  +  o»*+fl«— 0  (4) 

Based  on  the  Routh-Hurwitz  stability  condition,  the  real 
part  of  the  root  of  the  quadratic  polynominal  eigen  function 
(4)  is  negative  when  the  coefficients  a^,  a^,  a2,  a^  and  a^ 
i  have  the  same  sign  as  the  discriminant: 


/?a*d,aia*— a«o*— o«fl| 
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In  this  case,  the  motion  is  stable.  From  equation  (3)  we 


know  that  a^,  a^,  and  a^  are  always  greater  than 

zero.  Therefore,  R  should  also  be  greater  than  zero.  To 
this  end,  the  following  must  be  satisfied: 

>aia*  +  a4a*  (5) 


By  substituting  (3)  into  (5)  an  inequality  with  the  damping 
coefficient  h  as  ah  unknown  can  be  obtained.  This 
inequality  can  be  solved  with  respect  to  h.  This  is  the 
formula  needed  to  calculate  the  minimum  damping  coefficient. 


A»fc*+Aifc*+Atk+ A»=»o 


where 


Ai*a’»V+fljV* 

At-b'+brV'+btV* 

A,««iV+c,V* 

A,-d*V» 


(6) 


(7) 


In  equation  (7) 

a  a 

a\=ll<*lP{b+P-alXt-P-) 

a  a 

eOCM-af*) 


cl*sat(b+£-at)(fit  —  l)  +  (ati  +  b)')  .'  .. 

r,  -13)8* 

d.—n/J 

Now,  we  have  all  the  equations  required  to  solve  the  linear 

damping  coefficients  which  include  8  characteristic 

parameters,  where  I,  h  and  t  are  structural  parameters,  a, 

b, a  and  8  are  tire  parameters,  and  V  is  a  rolling  parameter. 

When  the  load  is  determined,  the  parameters  I,  t,  a,  b,  a  and  a 

are  also  fixed.  Hence,  other  parameters  can  be  calculated 

V 

from  equation  (8).  With  a  specified  speedy  the  parameters 
in  equation  (7)  can  also  be  determined.  Therefore,  all  the 
coefficients  in  equation  (6)  can  be  determined.  The  minimum 
damping  coefficient  h  can  be  obtained  by  solving  (6). 


direction  of  motionv 

\ 


center  of 
contact  plane 


center  of  direccional  axis 


wheel  axis 


Figure  1  Tire  Distortion  during  Wheel  Shimmy 
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II.  Assessment  of  Damping  Characteristics  of  the  Damper  and 
Treatment  of  Non>-linear  Problems 


The  minimum  damping  value  hmin  obtained  in  the  previous 
section  is  obtained  under  the  assumption  that  the  damper  is 
linear,  i.e.  the  damping  moment  is  proportional  to  the 
rotational  angular  speed.  In  reality,  most  hydraulic 
dampers  are  non-linear.  In  general,  the  relation  between 
the  damping  moment  M  and  the  angular  speed  of  the  driving 
lever  e  is 

M 0 1 6  * < ,  /g\ 


where  -  friction  moment,  which  is  produced  by  the 
friction  between  the  piston  and  the  inner  wall  of  the 
damper, 

a,  b  -  known  constants, 

0ti.  -  angular  speed  of  the  driving  lever  of  the 
damper. 

It  was  pointed  out  in  reference  [2]  that  the  drag  of  the 
fluid  is  determined  by  the  speed  at  which  it  flows  through  a 
channel.  The  drag  coefficient  of  the  fluid  through  a  hole 
is  determined  by  the  shape  of  the  jet  and  the  velocity  of 
the  fluid  in  the  jet.  If  the  channel  is  not  very  long  and 
the  speed  is  not  very  fast,  drag  is  proportional  to  the 
square  of  the  fluid  velocity.  If  the  channel  is  long  and 
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flow  speed  is  slow,  then  the  drag  is  proportional  to  the 

flow  speed.  The  damper  shown  in  Figure  2  belongs  to  the 

* 

former  type.  It  satisfies  the  following: 

In  the  M  -0*plot,  the  damping  coefficient  b  is  a  constant. 


Figure  2  Structural  Diagram  of  the  JXX  Damper 

Non-linear  dampers  have  the  following  advantages: 

(1)  It  is  easy  to  machine  the  flow  restrictor  hole  on  the 
piston.  It  will  not  be*  necessary  to  scrap  the  piston  when 
the  hole  is  drilled  too  deep  beyond  tolerance  or  the  drill 
bit  is  broken  inside  and  cannot  be  removed.  This  is  because 
this  kind  of  flow  restrictor  hole  is  shorter  and  the 
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machining  is  done  on  a  special  screw  plug  which  is  then 

screwed  onto  the  two  sides  of  the  piston  to  avoid  the 

* 

problems  described  above. 

(2)  The  non-linear  damper  is  more  effective  to  prevent 
shimmy.  It  increases  with  the  angular  velocity  of  the 
driving  lever.  The  damping  moment  picks  up  faster  than  that 
produced  by  a  linear  damper.  When  the  angular  velocity  is 
small/  the  damping  moment  is  smaller  than  that  generated  by 
a  linear  damper.  Therefore,  the  turning  performance  would 
not  be  affected.  Hence,  this  type  of  damper  is  recommended 
to  effectively  prevent  shimmy. 

However,  when  this  type  of  damper  is  used,  how  is  the  non¬ 
linear  damping  compared  to  the  theoretical  minimum  damping 
in  order  to  assess  its  ability  to  overcome  shimmy?  Two 
methods  are  introduced  in  this  paper.  Both  are  convenient 
and  reliable  to  use  in  various  occasions. 


»(!/**> 

<») 


(b)  damping  moment  vs. 

square  angular  velocity 
damping  aperture  diameter  1— 2—  M.s»n»  s— ♦i.amm  *— +>.omm 

5— ♦l.Smm  I— 4l.ta« 

Figure  3  Damping  Characteristics  of  the  JXX  Damper 


(a)  damping  moment  vs, 
angular  velocity 


1.  Minimum  Equivalent  Linearization  Method 


Let  the  yawing  angle  be  0  and  assume  it  is  in  simple 
harmonic  motion: 

0*Hsin©t 

(10) 


where  e  is  'he  angular  frequency  and  H  is  the  amplitude  of 
vibration. 


(1) 


energy  consumed  by  the  linear  damper  per  cycle 

r  2*  r*/2  , 

A—  |  Md0-4  (h  6  )(Ho>cosa>tdt) 

J  0  ■'o 


00 


9 


(2)  energy  consumed  by  the  non-linear  damper  per  cycle 

M«Ma+a*+fr0  * 

fin  e  n/l 

A'  =  Md0==4  *bd  *)d6 

J  o  J  0 


»4M,H+effo>H*+|-ioj*H*  (12) 

The  equivalent  damping  coefficient  is  determined  by  the  fact 
that  the  energy  consumption  per  cycle  is  equal.  Prom  A  *  A* 
we  get 


(13) 


Bence*  h  is  a  function  of  the  product  of  and  H.  Because 
we  are  looking  for  the  minimum  equivalent  damp  to  prevent 
shimmy*  we  have  to  find  its  extremum.  In  order  to  simplify 
the  process*  let  us  assume  p  =*»H  and  substitute  it  into 
113): 


(I3fl) 


From  9k/*>  =  0  we  get 


6=0 


Because  a  and  H  are  positive* 


P“ 


(M) 
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Also  because*  8th/8pt=!8M0/nptX),  h  has  a  minimum. 

Substitute  (14)  into  (13a) ,  we  get  the  formula  to  calculate 

► 

the  minimum  equivalent  linear  damping  coefficient 


h. 


i  • 


=  +  /3 Mj\ 

*  y  3M./26  4  3  V  2 i  / 


05) 


We  can  first  determine  the  constants  MQ,  a  and  b  by 
experimentally  measuring  the  damping  characteristics,  or 
performing  hydraulic  computation,  or  using  mathematical  fit 
and  then  find  the  equivalent  linear  damping  value  for  the 
non-linear  damper  using  (15). 


2 •  Tangent  Method 


This  method  does  not  require  to  find  the  constants  in 
equation  (9).  Instead,  the  data  obtained  in  test  the 
damping  characteristics  of  the  damper  are  plotted  on  the  M-d 
curve.  M  is  the  moment  applied  onto  the  damper  and  6  is  the 
angular  velocity  of  the  driving  lever.  A  tangent  is  drawn 
through  the  origin  o  and  the  slope  k  is  the  equivalent 
linear  damp  of  the  curve,  as  shown  in  Figure  4. 

The  damping  value  obtained  from  the  slope  k  of  the  tangent 
^  is  very  close  to  that  calculated  based  on  equation  (15) . 

The  significance  of  these  values  is  also  clear,  as  shown 
in  Figure  4. 
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+  «  i  +  *#* 


Figure  4  Equivalent  Linear  Damp  Obtained  by  Tangent  Method 

Let  the  equation  of  the  curve  and  the  equation  of  the 
tangent  passing  through  the  origin  be 


+b6* 

M-hd 


(16)' 


At  point  fit, 

M'+efi  +b6*=k6  .( 


i  .e. 


k**+(a-k)d  +M.-0  (17) 

* 

Let  us  solve  this  equation  and  we  get 

0  -~(fl~*)  ± 
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By  analyzing  (18)  we  know  that  only  when 


(a-k)*-4&M,=0  (19) 


and 


a— k<0  (20) 

the  straight  line  and  the  curve  can  intersect  each  other  in 
the  first  quadrant.  Moreover,  they  intersect  at  only  one 
point  (d„,  M„).  Proa  equations  (19)  and  (20)  we  get 

«-k— «/4Mfr 


Therefore, 


k«2  V^M.+a 


(21) 


Equation  (21)  is  an  expression  for  the  slope  of  the  tangent 
equation  in  Pigure  4.  The  coordinate  of  the  point  of 
tangency  is 


(22) 


In  order  to  explain  the  significance  of  the  slope  expressed 
by  (21),  we  calculated  the  constant  in  equation  (15). 


i.  _  8  /2 bM0  L  - 

hm"  *i-r-+a 


“2.0792v/frM#+a  (23) 

Comparing  (21)  to  (23),  we  arrived  at  almost  the  same  result 
using  two  different  approaches.  The  only  difference  is  in 
the  coefficient.  Therefore,  the  geometric  meaning  of  the 
non-linear  damper  derived  from  equation  (15)  is  equivalent 
to  the  tangent  going  through  the  origin.  Furthermore,  it  is 
unique • 

It  is  worthwhile  pointing  out  that  the  minimum  equivalent 

damping  value  obtained  with  either  method  is  very 

conservative.  However,  it  is  necessary  to  prevent  wheel 

shimmy  because  the  minimum  linear  damping  value  computed 

from  equation  (6)  is  the  critical  value  to  keep  the  motion 

stable.  In  practice,  we  must  consider  a  -t*i. 5^3.0  safety 

factor  to  ensure  a  safe  margin  for  stability.  Because  the 

(41 

upper  limit  that  would  affect  stability  is  very  high1  , 
this  minimum  method  is  certainly  feasible. 

III.  Example 

The  JXX  front  wheel  shimmy  damper is  used  as  an  example 
to  calculate  the  minimum  necessary  damping  value  to  satisfy 
the  entire  range  of  load  and  speed  based  on  equations  (6«*6). 
It  is  converted  to  kJJ*  *  14.1  Urns. 


14 


The  non-linear  damper  chosen  is  shown  in  Figure  2.  The  flow 

» 

restrictor  hole  was  determined  by  referring  to  that  in 
similar  airplanes  and  by  a  rough  estimate.  It  ranges  from 
approximately  j!2.0mm  in  diameter.  Figure 

5  shows  the  ft-0  curve  based  on  measured  damping  data. 


e 

x 


a  1  2  3  i  s 

9 

flow  restrictor  diameter:  1  -41.5  mm,  2  -41.8  mm,  3  -42.0  mm 
Figure  5  Damping  Characteristics  of  the  JXX  Damper 

Khich  damping  curve  do  we  choose  when  the  safety  factor  f  = 

1.8? 

He  first  used  the  tangent  method  to  find  the  solution.  The 
damp  value  is  designed  to  be 

fc;£. -1.8X14. 1  »25.4Nms 

The  tangent  lines  for  curves  1,  2  and  3  that  also  pass  the 
origin  are  ch^,  ch^  and  ch^.  Their  slopes  are  shown  in 
Table  1 
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Table  1  Slopes 


unit:  Nms 


M(N'm) 

0  (rad/s) 

K=M/  6 

curve  chj 

98 

2.6 

37.7 

curve  ch2 

98 

3.75 

26.1 

curve  ch^ 

98 

5.25 

18.7 

We  can  see  that  curve  ch^  meets  the  design  requirement.  The 
diameter  of  the  flow  restrictor  should  be<H.8  mm. 

Let  us  use  the  analytical  method  to  solve  the  problem.  The 

data  measured  with  the<t>1.8  mm  diameter  hole  is  used  to 

determine  the  constants  in  (9):  M  =  7.35  Nm,  a  =  0  and  b  = 

2 

23  Nms  .  The  result  obtained  from  equation  (21)  is 

*“2Xv'23X7.35*=26Nms 


When  using  the  minimum  equivalent  linear  damping  formula 
(15),  it  is 

The  relative  error  between  the  tangent  and  the  equivalent 
linearization  method  is 

27-2G 

c  » - 27 - *  100#  =-3.7# 

The  shimmy  experiments  and  long  term  use  of  the  JXX  landing 
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gear  have  demonstrated  that  this  assessment  method  is  safe 

and  reliable.  Moreover,  the  turning  ability  is  also  very 

* 

good. 
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Abstract 

This  paper  discusses  the  requirements  for  measuring  ground 
load  on  an  aircraft.  A  practical  data  treatment  method  has 
been  developed  for  the  case  that  the  measurement  satisfy  the 
minimum  requirement.  Although  the  method  is  illustrated 
with  taxiing  data  in  this  paper ,  however,  it  is  also 
applicable  to  the  treatment  of  turning  and  braking  load 
data.  Cht'n-e.  s  *  &  f  '  y  )  -rr— 

1.  Introduction 

/ 

Landing  gear  fatigue  is  a  serious  problem  to  consider. 

Among  the  Airworthiness  Directives  issued  by  the  FAA  between 
1955  and  1963,  41%  of  the  orders  for  piston  engine  airplanes 
i  and  31%  of  the  orders  for  jets  are  related  to  the  structure 
of  the  landing  gear. 1*1  More  accurate  estimation  of  fatigue 
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life  and  crack  propagation  will  require  a  more  accurate  load 
spectrum.  In  the  load  spectrum,  the  landing  and  takeoff 
taxiing  are  most  critical.1  *  1  Therefore,  the  focus  of 
discussion  is  centered  on  the  treatment  of  taxiing  load. 
Nevertheless,  a  practical  data  treatment  method  can  be 

r '  “ 

easily  used  to  treat  other  load  data. 

Statistical  treatment  of  the  taxiing  load  data  is  relatively 
complicated.  The  authors  did  a  preliminary  investigation  on 
this  subject  matter  in  reference  [6].  In  practice,  we  must 
consider  the  limitations  in  manpower,  financial  resources 
and  flight  measurement  conditions  and  take  a  series  of 
necessary  engineering  approximations.  For  example,  there 
might  be  some  correlations  between  the  various  parameters 
measured  and  they  are  very  important  in  the  force  analysis 
of  the  structure.  However,  they  are  difficult  to  obtain  and 
it  is  necessary  to  make  some  assumptions.^  On  the  basis 
of  reference  (6) ,  a  practical  statistical  treatment  of  the 
round  load  data  is  presented  in  this  paper. 

II.  Requirements  for  Measured  Results 

1.  Minimum  Requirements 

(1)  Overload  record  at  the  center  of  gravity  of  the 
aircraft  with  dual  parameter  counting,  and  the  frequency  at 
which  load  amplitude  and  mean  value  occur  at  the  same  time. 
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(2)  Loading  record  of  3  landing  gears  in  3  directions  with 

9 

dual  parameter  counting/  and  the  frequency  at  which  load 
amplitude  and  mean  value  occur  at  the  same  time. 

It  is  not  necessary  to  record  the  phase  correlation  between 
the  parameters  measured  above. 

2.  Satisfactory  Requirements 

(1)  The  comparative  relation  between  the  vertical  overload 
cumulative  frequency  curve  at  the  center  of  gravity  and  the 
vertical  load  cumulative  frequency  curve  of  the  landing 
gear.  Statistical  data  shows  that  the  primary  taxiing 
frequency  is.  3/4  Hz  and  could  be  as  high  as  2^4  Hz.  For 
small  aircraft  with  a  first  order  wing  frequency  of  7  ~  8  Hz 
it  is  considered  as  a  rigid  body.  The  comparison  mentioned 
above  could  provide  a  clear  correlation.  For  larger 
airplanes  this  correlation  does  not  exist  because  of 
flexibility.  A  direct  dynamic  analysis  is  required. 

(2)  In  the  landing  impact  stage/  the  correlation  between 
(front  and  back)  and  Py  (up  and  down)  are  recorded  and  a 
joint  distribution  is  made.  If  the  actual  data  is  not 
available*  use  the  statistical  data  shown  in  reference  [5). 
For  example*  the  landing  drag  coefficient  for  F-100  or  F-104G 
is  0.6*  the  rebound  coefficient  is  0.5  and  the  lateral 
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coefficient  (for  the  main  landing  gear)  is  0.5 


(3)  Record  the  relation  between  P  (left  and  right)  and  P 

z  y 

in  taxiing  and  make  a  joint  distribution.  If  there  is  no 

actual  data*  choose  lateral  overload  to  be  n  =  0.05  and 

z 

combine  with  n$  according  to  certain  frequency  and  phase  as 
suggested  in  reference  (5]. 


III.  A  Practical  Treatment  of  Taxiing  Data 

When  the  data  satisfies  the  minimum  requirements*  it  is 
possible  to  obtain  a  table  of  joint  occurrence  frequency 
involving  the  load  amplitudes  in  all  directions  and  the  mean 
values  at  the  center  of  gravity  or  on  the  landing  gear  (see 
Table  1) .  In  Table  l*  f  is  the  random  load  amplitude,  x 
represents  its  value*  <j  is  the  mean  load,  and  y  is  its 
value.  Z  and  *  are  in  unit  of  overload  g. 


Table  1  Joint  Occurrence  Frequency  of  Load  Amplitude  (f) 
and  Mean  Value  ( n  ) _ _ _ _ _ 

'N.  I  Ji  Ts  ~  Ti  —  ~ 
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The  joint  probability  distribution  of  a  two-dimensional 

random  variable  ( ]!,  rjj  )  can  be  estimated  based  on  the  data 

shown  in  Table  1.  On  this  basis,  the  variable  mean  load 
f  81 

spectrum1  J  can  then  be  derived.  However,  this  method 
requires  more  data  and  more  statistical  processes. 
Furthermore,  it  fails  to  emphasize  the  fact  that  load 
amplitude  is  a  major  factor.  In  reality,  from  the  load  data 
and  equivalent  life  curves  of  materials,  the  effect  of  load 
amplitude*  on  fatigue  life  is  far  more  greater  than  that  of 
mean  load  value.  The  following  method  takes  this  major 
factor  into  account  and  simplifies  the  treatment  for 
engineering  applications.  The  concept  is  to  first  convert 
all  the  mean  values  measured  to  normalized  mean  values  (such 
as  on  the  1  g  stress  level)  based  on  the  equivalent  life 
curve  and  then  treat  the  load  amplitude  statistically. 
Finally,  a  simplified  load  spectrum  with  a  given  order  is 
derived  from  the  normalized  damage  principle. 

The  purpose  of  measuring  the  ground  load  spectrum  of  an 
aircraft  is  to  provide  the  original  load  spectrum  to  the 
design  department.  This  is  a  standard  load  spectrum. 
Therefore,  we  cannot  use  the  fatigue  curve  (sjch  as  the  S-N 
curve  or  equivalent  life  curve)  of  a  specific  assembly  when 
we  analyze  the  measured  data.  Instead,  we  should  determine 
a  more  general  and  representative  fatigue  curve  and  then 
carry  out  the  necessary  simplification  procedures  based  on 
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the  measured  load  data. 

* 

1.  Several  Assumptions 

(1)  It  is  assumed  based  on  a  lot  of  data^  that  the  stress 
amplitude  S_  and  fatigue  life  N  have  the  following  relation 

when  the  mean  stress  S_  is  a  constant: 

■ 

IgN^c+MgS.  (  i  ) 

where  a  and  b  are  material  constants. 

(2)  The  equivalent  life  S_  -  Sm  curve  is  a  straight  line 

d  m 

which  passes  through  S_j)  and  (Sb»  0).  The  equation  is 

s-“s"<'- ir>  <J) 

where  Sb  is  the  stretching  limit  of  the  chose  material  and 
is  the  endurance  limit  of  the  material  under  a  symmetric 
cycle  (R=-I ) . 

(3)  The  fatigue  damage  follows  the  Miner  linear  cumulative 
damage  law. 

(4)  When  the  design  load  coefficient  n  is  reached,  the 
total  stress  at  the  ideal  endangered  point  reaches  the 
stretching  limit  S^.  Based  on  this,  the  stress  level  at  1  g 


23 


during  taxiing  is 


Si. 


„  S, 

X*»ts 


2.  Conversion  to  Equivalent  Life 


(3) 


The  data  in  Table  1  are  converted  to  dual  parameter  stress 
at  the  endangered  point. 


— Si* *X(,  i”l ,  2*  *** >  ^  | 

ml  (<) 


The  dual  parameter  stress  (S  .  ,  S  . )  has  the  same  frequency 

a  X  m  j 

of  occurrence  as  that  of  ( x yj),  ni j • 


Based  on  the  equivalent  life  principle,  the  dual  parameter 

stress  data  are  simplified  to  single  parameter  stress  data 

with  a  constant  mean  stress  value.  The  constant  mean  stress 

is  denoted  as  Sffi0.  If  the  equivalent  life  conversion  of  any 

pair  of  stress  (S,.  ,  S  . )  is  (S  ..  ,  S^) ,  then  based  on 

a  l  n  j  a  i  3  mo 

assumption  (2)  we  have 


i“l 

1*1  pin 


(5) 


The  frequency  of  occurrence  of  S 


ls  stm  nij 
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3.  Equivalent  Damage  Conversion 

* 

After  equivalent  life  conversion,  the  order  of  stress 
amplitude  is  usually  still  too  high  to  use.  Therefore,  we 
still  have  to  make  another  conversion  based  on  the 
equivalent  damage  principle. 

(1)  Direct  Conversion 

The  order  of  converted  stress  amplitude  is  pre-set  at  M 
based  on  pre-determined  requirements.  Choose  M  -  8  for 
taxiing  load  by  referring  to  reference  [3]  Denote  that  the 
stress  amplitude  at  each  order  is  S'i*  (k  =  1,  2,  ..,  M). 

Sai^  is  converted  to  its  closest  stress  amplitude  s^1. •  Let 

us  assume  that  the  conversion  order  is  nil’ .and  k^  is  the 

sequence  of  Sa^  to  be  converted  to  the  kth  order  of  stress 

amplitude.  N. .  represents  the  life  of  S  . .  at  the  ideal 
i  j  ai  j 

(k) 

endangered  point  and  B  is  the  life  of  Si*’  at  the  ideal 
critical  point.  Based  on  assumption  (1)  we  have 


N«|— !0*(S.»j)* 


i 

Again*  based  on  assumption  (3) 


Nrrr  N]t 


(8) 

(60 


(7) 


25 


From  aquations  (6).  (6*)  and  (7)»  it  is  not  difficult  to  get 


(8) 


T  * 

Add  all  the  numbers  simplified  to  the  same  order  of  stress 
Si«  and  we  have  the  number  of  times  n  that  stress  occurs 
v  in  a  landing  or  takeoff,  i.e. 

*»-» 


c< » >  j.c<  »  ) 

where  m^  fa  the  number  of  Sai..  occurring  in  (0,  — — —  ). 


•  2  - ) 

and  h  *  1,  ••••>,  B— !• 


•«*»-  2 

where  a^  fa  the  number  of  s»»j  that  is  greater  than  — 2 — jj - 


A  stress  amplitude  spectrum  with  a  constant  mean  value  S 

mo 

is  obtained  (see  Table  2). 


The  overload  coefficient  spectrum  with  a  mean  overload  of 

ST  /ST  can  ha  obtained  by  dividing  stress  with  S.  . 
mo  ig  lg 
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Table  2  Stress  Amplitude  Spectrum  with  a  Constant  Mean 


Value  Saij 


stress  amplitude  s«l> 

js«*>j~|s<-> 

- 

S‘*> 

.  *  *>!... ! 

number  of  occurrences  I  (  1 

|  a*  *  *|  •“  l  n<u> 

(2)  Statistical  Smoothing  Conversion  Method 

The  equivalent  life  converted  stress  amplitude  can  be 
smoothed  out  by  using  an  appropriate  probability 
distribution  function  before  performing  the  equivalent 
damage  conversion. 

(a)  selection  of  suitable  stress  amplitude  probability 
function 


In  order  to  select  a  suitable  function  from  several 
probability  distributions*  an  optimization  check  for  the  fit 
can  be  used.  To  this  end*  the  Sai^  obtained  are  rearranged 
from  low  to  high  (as  shown  in  Table  3)*  where 


ft' 

2*» !>**«* 
l»l 


(10) 


27 


Table  3  Table  of  Sajj  in  Sequence 


stress 

amplitude 

i  1  |  1  1  1 

|  S«<  »>j  S«c*>  ~j  S«<l>  i~j  S«  <•'  ) 

number 

of  occurrences 

!■<»>  )*<*>  j  a < at  > 

There  are  many  methods  to  check  the  fit.  Here,  we  are  using 
the  W  method  introduced  by  Simlnov: 

Ol) 


where  F.  and  F  are  the  empirical  and  theoretical 
no 

distribution  functions  of  Sa4  respectively. 

2 

n*W  can  be  calculated  based  on  the  following  equation: 


•'  »c»> 

**Hr*~-rfer+2  2c«.  02) 

***  im\  hm\ 

where 

[f.(S\u)~-^^— j  »*!, 

Gu*  ,*  **  > 

z(S*o>)+2*-I 

[fiCS'ii*)  2n*  “]»  *^1. 

o  ii(  i>  =0 
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3*c  *>_ 


■S«<i i 


2n<  i 


i  =  l 


~frS«u-i >  -|.f.3a<ut)~S«<<-i>  j <i<»' 

2  2n « i > 


}  (U) 


S«/m-  >  -  i  )  .fr.fr.  S«<«.  >  S«<»'  -  i  >  j  =  n' 

«(•’  > 


2 

Based  on  the  given  degree  of  significance  a ,  ,  the  percentage 

,2 


2.  corresponding  to  100  a  t  can  be  looked  up  from  the  n*W4 

percentage  table.  If  n*W  <  zt,  ,  FQ  is  considered  as  the 

probability  distribution  function  of  S  . ..  If  several 

ai  j 

theoretical  functions  pass  the  test,  the  one  with  the 
2 

smallest  n*W  can  usually  be  chosen  as  the  distribution 

function  for  S  .  .  • 
aij 


Discrete  Stress  Amplitude 


The  range  of  Sai..,  {&Q,  dk2>,  is  divided  into  *2  intervals 

of  equal  width.  Based  on  the  probability  function  Fo(x)  of 

S  . •*  we  can  calculate  the  probability  of  occurrence  of  S  . . 
ax  j  sxj 

in  each  interval  and  use  it  as  the  probability  for  the 
stress  amplitude  to  occur  at  the  center  of  that  interval. 
Based  on  the  number  of  loadings  in  one  takeoff  or  landing, 
the  frequency  for  the  stress  amplitude  to  occur  at  the 
center  of  each  interval  can  be  calculated.  The  width  of 
each  interval  is 

W-iisf*!  (is) 

In  the  ith  interval  (dj_j.,  di)  (i  *  1,  2,  ...,  k^),  the 
representative  stress  is  denoted  as  S  . : 
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(16) 


S.,=dt  +  (i-y)W 

The  probability  for  S  •  to  occur  is 

u  X 

p,— F.(di)—F*(di-i) 

In  one  landing  and  takeoff,  the  number  of  occurrence  of  S  , 

d  x 

is 

**“»*ln  (is) 


(c)  The  method  used  to  perform  equivalent  damage  conversion 

of  the  statistically  treated  data  is  similar  to  that 

described  before  in  (1).  S_..  and  n . .  are  converted  to  S  . 

aij  ij  ai 

and  n^#  respectively  to  get  the  simplified  stress  spectrum 
and  overload  coefficient  spectrum. 

No  matter  which  method  we  use,  after  the  simplified  stress 
spectrum  and  overload  coefficient  spectrum  are  obtained,  we 
can  derive  the  cumulative  frequency  spectrum  for  stress  or 
overload  coefficient  based  on  need. 

4.  Statistical  Treatment  of  Data  on  Multiple  Landings  and 
Takeoffs 

The  taxiing  load  data  were  obtained  in  L  landings  and 
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takeoffs  under  nominally  identical  conditions.  The  taxiing 
data,  for  each  run  can  be  processed  to  obtain  the  simplified 
overload  coefficient  spectrum.  The  results  are  shown  in 
Table  4. 


Table  4  Simplified  Overload  Coefficient  Frequency  Table  for 
Multiple  Landings  and  Takeoffs 


/V) 

CD  If  for  an  overload  coefficient  X  (k  *  1,  2, —  ,  M), 
the  frequency  of  occurrence  in  each  taxiing  follows  the 
normal  distribution  (which  was  found  to  be  valid 
experimentally)#  it  is  possible  to  calculate  the  upper  limit 
of  confidence  for  the  mean  occurrences  of  an  overload 
coefficient 

c<  »  > 

-7== 

09) 
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where 


*“>=*  rl-2«i‘> 

**  i  •  1 

i»il  J 


t  is  the  point  corresponding  to  100  a  %  in  the  t 
distribution  which  has  v  =  L-l  degrees  of  freedom.  It  is 
looked  up  from  the  i-a distribution  table  based  on  the  given 
confidence  level. 


Finally,  the  upper  limit  spectrum  for  each  overload 
coefficient  in  a  single  taxiing  event  is  obtained  (see  Table 


Table  S  Upper  Limit  of  Mean  Humber  of  Occurrences  for  Each 
Overload  Coefficient 


overload  coefficient 


number  of  occurrences  I a* 


(2)  Based  on  need,  it  is  also  possible  to  calculate  the 
upper  confidence  limit  for  the  number  of  occurrences  for 
each  X^  in  each  landing  or  takeoff,  . 


where  n^l,  and  t.  are  defined  as  before. 


Thus#  the  upper  limit  spectrum  for  the  overload  coefficients 
to  occur  in  one  taxiing  is  obtained. 

* 

With  the  above  result#  the  cumulative  frequency  spectrum  of 
the  overload  coefficients  in  one  taxiing  can  be  conveniently 
derived. 

It  is  easy  to  see  that  the  above  load  treatment  method  can 
be  easily  expanded  to  treat  turning  and  braking  data 
measured  on  the  ground. 

IV.  Discussion  and  Conclusions 

In  addition  to  the  equivalent  damage  method#  which  was 

described  earlier  to  convert  the  equivalent  life  data  (i.e. 

number  of  occurrences  for  each  stress  amplitude  of  a 

constant  mean  value)  to  the  effective  number  of  occurrences 

of  each  simplified  stress  amplitude#  there  is  another 

equivalent  damage  probability  method  which  guarantees  equal 

probability  of  stress  induced  damage  before  and  after 

f81 

conversion.  However#  it  can  be  verified  that  when 
assumption  (1)  is  valid  the  fatigue  life  follows  a 
logarithmic  normal  distribution  at  all  stress  levels.  The 
standard  deviation  is  a  constant.  (It  was  proven 
experimentally  that  aeronautical  metals  can  usually  meet 
this  requirement.)  In  this  case#  the  equivalent  damage 
probability  method  is  essentially  the  same  as  the  equivalent 


33 


damage  method.  In  other  words,  the  equivalent  damage  method 

used  in  this  work  is  the  equivalent  damage  probability 

» 

method  under  less  rigorous  conditions. 

The  treatment  of  aircraft  ground  load  data  is  a  very  complex 
job.  The  method  presented  in  this  paper  is  a  statistical 
treatment  of  the  data  under  the  premise  that  the  minimum 
requirements  are  met.  Figure  1  shows  the  flowchart  of  this 
method.  If  the  data  meets  the  satisfactory  requirements  in 
the  future,  new  statistical  methods  can  be  developed  based 
on  the  correlation  of  the  parameters.  « 

Statistical  treatment  of.  ground  load  data  is  only  a  portion 
of  the  study  of  ground  load  spectrum  and  the  formulation  of 
standards.  The  result  does  not  represent  a  standardized 
load  spectrum.  Many  factors  such  as  load  conditions  and 
environmental  impact  are  difficult  to  be  reflected  in  finite 
measurements.  They  must  be  considered  by  other  suitable 
methods  (such  as  empirical  parameters). 

The  methods  introduced  in  this  paper  has  been  successfully 
used  to  treat  the  data  obtained  in  120  landings  and  takeoffs 
of  the  Chinese  made  Yun-7  transport  at  different  airports. 

A  landing  load  cumulative  frequency  spectrum  similar  to  the 
one  described  in  MIL-A-008866B(USAF)  has  been  obtained.  In 
addition,  it  is  being  used  in  the  fatigue  testing  of  the 
landing  gear  of  that  aircraft. 


Figure  i  Flowchart  for  the  Statistical  Treatment  of  Data 

1.  ig  stress  level,  S(  a ,  at  representative  danger  pciiU  i 
the  cross-section;  2.  family  of  equivalent  life  curves; 
3.  measured  load  data;  4.  dual  parameter  overload 
coefficient  spectrum;  5.  dual  paj  afneter  strews  spectrum 
representative  danger  point  in  the  cross- sec  Li  on  ;  <:•. 

equivalent  life  conversion;  7.  constant  mean  stress 
amplitude  spectrum;  3.  S-N  curve  cumulative  damage-  tiiii; 
9.  direct  equivalent  damage  conversion;  13.  simplified 
stress  amplitude  spectrum;  11.  simplified  overload 
coefficient  spectrum;  12.  upper  limit  frequency  spe^frui., 
for  the  occurrences  of  simplified  overload  coefficients; 
13.  fitting  stress  amplitude  with  a  distribution  function 
*4.  equivalent  damage  conversion;  15.  simplified  stress 
amplitude  spectrum;  16.  simplified  overload  coefficient 
spectrum;  17.  upper  limit  frequency  spectrum  for  the 
uccurrem-ea  of  simplified  overload  coefficients;  18. 
equivalent  damage  probability  conversion;  IS.  simplified 
stress  amplitude  spectrum;  20.  simplified  overload 
coefficient  spectrum;  21.  upper  limit  frequency  spectrum 
for  tire  occurrences  of  simplified  overload  coefficients; 
22.  S-N  curve  fatigue  life  probability  distribution 
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